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2 ■ We determine the intrinsic shapes and orientations of 27, 450 type I and II active 

galactic nucleus (AGN) galaxies in the spectroscopic sample of the Sloan Digital Sky 
Survey Data Release 7, by studying the distribution of projected axis ratios of AGN 
hosts. Our aim is to study possible alignments between the AGN and host galaxy 
systems (e.g. the accretion disc and the galaxy angular momentum) and the effect of 
dust obscuration geometry on the AGN type. We define control samples of non-AGN 
galaxies that mimic the morphology, colour, luminosity and concentration distribu- 
tions of the AGN population, taking into account the effects of dust extinction and 
, reddening. Assuming that AGN galaxies have the same underlying three-dimensional 

■ shape distribution as their corresponding control samples, we find that the spiral and 

elliptical type I AGN populations are strongly skewed toward face-on galaxies, while 
ellipticals and spirals type II AGN are skewed toward edge-on orientations. These find- 
ings rule out random orientations for AGN hosts at high confidence for type I spirals 
(Sx 2 ~ 230) and type II ellipticals (<5% 2 ~ 15), while the signal for type I ellipticals 
and type II spirals is weaker (<5% 2 ~ 3 and S\ 2 ~ 6, respectively). We obtain a much 
stronger tendency for the type II spirals to be edge-on when just high [OIII] equivalent 
K> , width (EW) AGN are considered, suggesting that > 20% of low [OIII] EW edge-on 

j_j ■ type II AGN may be missing from the optical sample. Galactic dust absorption of the 

broad-line region alone cannot explain the observed inclination angle and projected 
axis ratio distributions of type I and II Seyfert types, implying that obscuration by a 
small-scale circumnuclear torus is necessary. These results favour a scenario in which 
the angular momentum of the material which feeds the black hole retains a memory 
of its original gas source at least to some small, non-negligible degree. 
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1 INTRODUCTION 

Galaxies with Active Galactic Nuclei (AGN) have become 
a major focus in extragalactic studies due to their role in 
galaxy formation scenarios favoured today. Feedback from 
the AGN, and its effects on gas infall and star formation are 
invoked to explain observed trends on galaxy st ar formation 



rates , luminosities and colours (see for instance . I Bower et al 



20061 : ICroton et al.ll2006l : jCattaneo et alJl20od:ISiiacki et al 

200' 



Hopkins et~al1 120071; iMarulli et al.l 120081 ; lLagos et al 



20081 : ISomerville et al.ll2008T l. However, many aspects of the 



relationship between AGN and their hosts remain unclear, 
such as the relation between the angular momentum of 
the AGN system (i.e. accretion disc and black hole spin) 
and the galaxy kinematics. A level of coherence between 
the two co uld influence the d evelopment of the black hole 
(BH) spin (|Lagos et alj|2009l ) which regulates the mass-to- 
ener gy conversion in radiatively efficient accretion p hases 
(e.g. Ikawlings fc Saunder3ll99ll : iMarconi et aTll2004r ), and 
which has be en postulated as ke y to the process of jet pro- 
duction (e.g. ISikora et al.l 12007^ : therefore, this could di- 
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rectly affect the energetic capacity of a BH to produce ef- 
fective feedback. 

A way to study galaxy- AGN orientations, for example in 
weak objects suc h as radio-quiet AGN (i.e. Seyfert galaxies, 
lOsterbrockll 19841 ). is via the measurement of the inclination 
angles relative to the line of sight of AGN host galaxies, 
since narrow (Seyfert II) and broad line (Seyfert I) AGN 
are thought to be the same physical phenom enon seen un- 
der different orientations fsee lAntonuccilll993l for a review) . 
This is called the " Unified AGN model" (e.g. lBarthellll989l; 
Madau et al. I ll994l ; lUrrv fc Padovanill 19951 ; |G unn & Shanks 



19991 ) in which the emission lines come from the AGN sys- 



tem composed of a BH surrounded by an accretion disc and 
an obscuring torus. The accretion disc is responsible for the 
ionizing spectrum of the AGN, and the broad line region 
(BLR) gas is distributed on similar scales; the direct ob- 
servation of these regions results in a Seyfert I galaxy. The 
opaque torus absorbs the optical emission from the accretion 
disc and the BLR along certain lines of sight, giving rise to 
a Seyfert II galaxy. 

The main evidence in favour of the unified model 
comes from observations of broad emission lines in polarised 
light from Seyfert II galaxies (e.g. [Antonucci fc Millet 



Mill er et al.lll990l; IZakamska et al.ll2005l ; IZhang et al 



Borguet et aIT l2008). and is further supported by e.g. 
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2008 



surem ents of similar BH accretion rates i n the two AGN 
types (|Netzerll2009l ; I Hatziminaoglou et al.l I2009T ) . However, 
the unified model does not imply a relation between the in- 
clination of the host galaxy and that of the obscuring torus 
(which may be linked with the BH sp in; IVolonteri et al.l 
120071 . iKing et al.ll2008l ; lLagos et al.ll2009h . leaving open the 
question of the existence of such a relation. Indications of a 
relation between the kinematic properties of the outer and 
inner galactic co mponents have been recently found in the 
S AURON survey jBacon et all l200ll ; Ide Zeeuw et all |2002| ) 
bv lDumas et al.l (|2007 ) in AGN Seyfert galaxies. These show 
differences in the position angle of outer and inner compo- 
nents typically of w 30°, thus these two components are not 
randomly oriented with respect to one other. In addition, 
type II AGN spirals appear to be more el ongated in optica l 
bands than are non-active spiral galaxies (|Shen et al.ll201oT ) , 
possibly indicating a tendency towards edge-on orientations 
in this type of AGN galaxies, which within the unified AGN 
model would point to a connection with the to rus orienta- 
tion. On the other hand, iKrai novic et al. I (|2008l ) found kine- 
matically decoupled cores (i.e. central regions) in 30% of 
a small sample of elliptical galaxies observed as part of the 
SAURON project. Thus, the evidence for alignment between 
the galaxy as a whole and the AGN is not conclusive. 

Most of the studies in the literature on the orientations 
of AGN have been done using powerful radio objects, by 
studying the level of alignment of relativistic jets and their 
hosts. Since the jet direction is set by the angular momen- 
tum of the accreti on disc (e.g. Lyndcn-Bcll 2006) and/or 
the BH spin (e.g. iBlandford fc Znaiekl 1 19771 ; iFragile et al.l 
120071 ; [b arausse fc Rezzollall2009l ). it can be directly com- 
pare d with the orientation of the host galaxy. Several papers 
(e.g.|Kinnev et al.|[2000l ; ISchmitt et~ail2002i ; iGreenhill et all 



120091 ) have found that radio and optical position angles 
are uncorrelated, although their samples contain < 100 ob- 
jects, meaning tha t the statistical precision of these studies 
is limited. Indeed, iBattve fc Browne! j2009) found a signif- 



icant level of alignment between the minor axes of a sam- 
ple of 6, 053 elliptical galaxies and the position angles of 
relativistic jets emitted from their central AGN. However, 
such an apparent alignment could be difficult to distin- 
guish from the "alignment effect" jMcCarthv et alj 1 19871 : 
IChambers et al.ll 19871 ; iMcNamara et al.lll996T and 1996b) in 
brigh t radio sources, in which the jet trigg ers star formation 
(e.g. |Pe Youndll995l ; iBlundell etahll 19991 ) or where the blu e 
component is simply scattered AGN light (|McCarthvlll993l ). 

The BH spin development is intimately linked with the 
growth of the BH itself. This means that its value and its 
level of alignment with the galactic components give im- 
portant cl ues about the way in which the accretion pro- 
ceeds (e .g. Moderski et al.lll99 1; iHughes fc Blandford| [2003; 



IShapird E005F The spin of a BH depends on whether it 
gained most o f its mass via mergers with other BHs or via 
accretion (e.g. Volonteri et al., 2005; Be rti fc Volonterlll2008l ; 
lLagos et al1l2009l ). IVolonteri et al.f(|2005l ) found that gas ac- 
cretion affects the spin evolution more than do mergers be- 
tween BHs, since binary coalescences alone do not lead to 
a systematic spin-up or -down in time, while a prolonged 
accretion period efficientl y spins BHs up t o the maximum 
value. On the other hand, IKing fc Pringkj |2006l ) suggested 
that the accretion proceeds via a series of short episodes as 
small amounts of gas with randoml y oriented angu l ar mo - 
menta fall in. In contrast with the IVolonteri et alj (|2005l ) 
model, in this scenario BHs can achieve only low spin val- 
ues, a= cJbh/G JV?bh ^ 0.1 — 0.3 through gas accretion 
(King et al. 2008) a nd r elatively high spin values (i.e. a< 0.7; 
iBerti fc Volonteril [2008) can only b e achieved through th e 
effects of galaxy and BH mergers (|Fanidakis et al.l l2009h . 
X-ray observations reveal relativistically broadened Fe-Ka 
fluorescence lines in several AGN, indicative of a> 0.9 



dlwasawa et al.lll996l;lFabian et aLll2002l ; lRevnolds fc Nowakl 
2003; Brenneman fc Reynolds! 120061 ). Another approach to 
constrain BH spins is via quasar demographics. The ob- 
served light from quasars is directly proportional to the 
rate at which supermassive black holes are accreting ma- 
terial, thus a comparison of the quasar luminosity function 
with the prese nt-day BH ma ss function constrains the radia- 
tive efficiency (|Soltanlll982l ). which in turn depends on BH 
spin. Analys es based on t his ap p roach have come t o different 
conclu sions: lElvis et all (|2002l) , ICao fc Lil (|2008l ). lYu fc Lul 
(2008) reported large energy conv ersion efficiencies (> 
13%), indicative of large spins, while lYu fc Tremainell2002l . 
iMarconi et~aT]|2004 iMartfnez-Sansigre fc Tavlorl (|2009l ). ~re- 
ported efficiencies as low as 7%, indicative of slowly rotating 
BHs. This may depend on the redshift, given that the typi- 
cal accretion rates of relatively ma ssive BHs are likel y to be 
much larger at earlier epochs (e.g. IWang et al. I l2009l ). 

Lagos, Pa dilla fc Cora (2009, LPC09 hereafter) used the 
a-disk model (IShakura fc Sunvaevl 1 197*31 ) to follow the BH 
spin devel opment in the sem i-analytic model of galaxy for- 
mation of lLagos et a l. (2008) within a cosmological frame- 
work and showed that, in the hypothetical case in which 
accretion discs preserve the angular momentum direction of 
the bulk of the cold gas involved in star formation, massive 
galaxies should host BHs with spin values close to unity, re- 
gardles s of additional ph ysical effects such as accretion disc 
warp s (|King et al.| [2005) or fragmentation (|King fc Pringld 
2006). We refer to this scenario, in which the BH is spun 
up smoothly via steady accretion of gas that conserves the 
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angular momentum direction of the original gas source, as 
the coherent model. Using the same semi-analytic model, 
the authors tested the effect of assigning random orienta- 
tions to the incoming material and show that high BH spin 
values can only be obtained by aligning the inner regions 
of the accretion disk with the BH spin (i.e. disc warps; 
iKing fc Pringkll2006l) ; this case will be referred to as the 
chaotic model. 

In this work we focus on the chaotic/coherent accretion 
dichotomy, by studying intrinsic shapes and inclinations of 
host galaxies of AGN with obscured a nd unobscured B LR in 
the Sloan Digital Sky Survey (SDSS.lYork et al.ll200oT ) Data 
Release 7 (DR7, lAbazaiian et al.1 12009| ). to determine the 
degree of alignment between galaxy and accretion discs (or 
tori). Throughout the paper the term AGN will refer only 
to Seyfert galaxies. We will assume that the torus and the 
accretion disc are aligned, with the exception of a possible 
warp in the a ccretion disc in th e vicinity of the central super- 
massive BH llKing et al. l l2005T ) as seen in AGN masers (e.g. 
IZhang et al.ll200a ). We compare the results obtained with 
predictions of the LPC09 model in the coherent and chaotic 
scenarios to constrain the physics of gas inflow. 

This paper is organised as follows. In Section 2, we 
describe the sample and the selection of control samples 
matched to the properties of AGN galaxies. We charac- 
terise the three-dimensional i ntrins ic shapes of these samples 
following IPadilla fc Strauss! ((2008). In Section 3 we study 
the orientations of AGN host galaxies, systematic effects in 
the AGN selection, and the origin of biases toward face- 
or edge-on orientations by analysing the effect of galactic 
disc and nuclear torus obscuration. In Section 4 we compare 
the observed inclination angle distributions with theoretical 
predictions for the chaotic and coherent models of LPC09. 
We discuss and summarise our main results and their im- 
plications for AGN models in Section 5 and Section 6, re- 
spectively. Throughout this paper we assume a cosmological 
model characterised by matter and dark-energy density pa- 
rameters Q m ~ 0.27, Qa = 0.73 and a Hubble constant 
Ho = 73kms~ Mpc -1 . These parameters are consistent 
with the results from the Wilkinson Microwave Anisotropy 
probe fWMAP. iHinshaw et al.1 120091 . iDunklev et~aH |2009| . 
ISanchez et alj|2006h . 



2 THE INTRINSIC SHAPES OF AGN HOST 
GALAXIES IN THE SDSS 

In this Section we characterise the three-dimensional intrin- 
sic shapes of the hosts of AGN. Our AGN s ample was iden- 
tified using the method of lHao et all (|2005ah from the SDSS 
DR7 spectra. We assume that the AGN have the same in- 
trinsic three-dimensional shapes as well-defined control sam- 
ples of non-AGN galaxies. Galaxies in the control samples 
should be free of inclination biases, and therefore their ob- 
served projected shapes will only reflect their intrinsic three- 
dimensional shapes, plus the effect of extinction. Compari- 
son of the axis ratio distribution of AGN galaxies with those 
of the control sample allows us to make inferences on the ori- 
entation distribution of the AGN. 



2.1 The SDSS DR7 AGN sample 

The AGN data used in this work consists of 6, 153 type I and 
21, 297 type II AGN identified by their emission-line proper- 
ties from about 698 , 000 galaxy spectra (the MAIN survey; 
IStrauss et ail |2002h in the SDSS DR7 using the methods 
of lHao et al.l (l2005al ). In order to exclude objects with low 
S/N emission lines, Hao et al. reject from the sample galax- 
ies whose Ha rest-fra me equivalent width (EW ) is < 3A. We 
exclude type III AGN (|Kauffmann et al.ll2003l), which lie be- 
twe en the star-forming and AGN loci ( Kewlev et al.ll200lT ) in 
the IVeilleux fc Osterbrockl (|l987f ) line ratio diagram, since 
much of their emission could come from star-formation ac- 
tivity. Hao et al. fit for and subtract the stellar continuum 
before measuring the strength and widths of emission lines. 
Ha and H/3 are fit to two Gaussians; those objects with 
a significant broad component (FWHM > 1200km/s) are 
termed type I, while thos e objects lying beyo nd the star- 
forming line ratio limits of iKewlev et al.l (|200ll ) are termed 
type II. 

We separate the AGNs further by morphology, using 
the fraction of r-band light which fits a de Vaucouleurs pro- 
file, ,fdev (from model fits available from the SDSS pipeline, 
lAbazaiian et al.l 120041 ). Early-type galaxies are selected by 
requiring /dev ^ 0.9, and we classify all other galaxies as 
spirals. We therefore define four AGN samples, selected by 
AGN type and morphological host galaxy type. 

We apply a dust extinction and reddening correction 
to luminosities and colours, respectively, of spirals using the 
model by Padilla & Strauss (2008, PS08 hereafter), where it 
is assumed that the amount of reddening is proportional 
to the path length of the light through the galaxy. The 
extinction is a function of the inclination angle 6 for any 
given galaxy, the mean galaxy height to diameter ratio, 
7=1 — C/B (where galaxies are modelled as triaxial el- 
lipsoids of major, middle and minor axis A, B, and C, re- 
spectively), and the edge-on extinction, Eo (in magnitudes). 
Here, a perfect edge-on galaxy will have 6 — 90°, while a 
face-on galaxy has 9 = 0°. For each spiral galaxy we ob- 
tain _Eo and 7 from Table 4 of PS08 (interpolating in abso- 
lute magnitude) . In turn, 6 is determined from the observed 
projected axis ratio by interpolating in their Table 8. The 
edge-on reddening in r magnitudes, _Ro, is related to Eo via 
-Eo = f Ro, where / = 2.77 is the reddening parameter. PS08 
showed that the model is not strongly dependent on /, al- 
lowing us to use this simple approximation to correct both 
colours and luminosities. The extinction affects the maxi- 
mum volume, Knax, out to which a galaxy can be detected 
in the SDSS flux-limited catalogue. We select galaxies with 
M r — 51ogio h < —17 in order to avoid very small values of 
Knax, since fainter galaxies tend to dominate the noise in the 
estimate of the distribution functions we consider below. 



2.2 Selection of Control Samples 

We select control samples matching the colours (g — r), lu- 
minosities (M r ) corrected for reddening and extinction, re- 
spectively, and concentration^] of the AGN host samples. 



1 Concentration is d efined as the ratio between the radii enclosing 
90% and 50% of the iPetrosianl lll976h light, r 9 o/r 5 o. 
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Table 1. Number of galaxies in each AGN sample and their cor- 
responding control sample. 



Sample 


AGN sample 


control sample 


Type I ellipticals 


2,712 


50,537 


Type II ellipticals 


7,653 


45,253 


Type I spirals 


3,360 


57,010 


Type II spirals 


13,578 


66,626 



All magnitudes shown are fc-corrected to z = 0.1. We then 
construct a grid in the parameter space denned by these 
three galaxy properties, and populate each cell with the 
AGN hosts weighted by 1/Knax- The control samples are 
constructed by filling each grid cell with randomly chosen 
main sample galaxies from the SDSS so that E V~^ x Gals w 
NxYl V^~ ax AGN ; N is chosen so as to maximise the size of the 
control samples, with a maximum departure from the target 
distribution in any grid cell of less than 20%. N is 5 — 20 for 
our four AGN samples. Notice that this selection takes into 
account any intrinsic correlation between these quantities in 
the AGN sample. This ensures that control samples mimic 
the AGN hosts as closely as possible. TableQ]shows the num- 
ber of galaxies in each AGN sample and their corresponding 
control sample. 

Fig. [T] shows the rgo/rso (left panel), M T — 51ogio h 
(middle-left panel) and g — r (middle-right panel) Knax- 
weighted distributions of each of the AGN subsamples 
(solid lines) and their corresponding control samples (dashed 
lines). The right panels show the redshift distributions of 
both populations. The control sample d istribution s all fol - 
low the AGN distributions very closely. lHao et all |2005b) 
noted that sufficiently low luminosity AGN would not be 
selected spectroscopically due to (i) undetectable emission 
lines, (ii) higher noise in the measurements of emission line 
ratios for weaker emission line strengths, or (iii) total galaxy 
fluxes (including the nuclear emission) below the spectro- 
scopic sample limit. This could result in a difference in the 
redshift distribution of AGN and non-AGN galaxies, which 
would cause a difference in the values of V ma x- The observed 
agreement between the redshift distributions of AGN and 
non-AGN galaxies (right panel of Fig. [1} suggests that this 
is not a problem. We also find a broad agreement between 
the M r — 51ogio h distributions of AGN and control samples 
for the spiral population, and only a small offset at the flux 
limit (m r « 17.5) in the elliptical population. We conclude 
that our V m &x does not need extra corrections for these ef- 
fects and is adequate for our analysis. 

In Fig. [TJ the host galaxy properties of the two AGN 
types (for the same morphology) show very similar photo- 
metric properties, a result consistent with the AGN unified 
model, in which the difference between the types is due solely 
to orientation. Even though the non-obscured nuclear emis- 
sion from type I AGN could affect the photometry of the 
galaxies, the similar distributions of luminosities and con- 
centrations of type I and type II AGN indicate that such 
an effect is minor. In fact, type I AGN galaxies are only 
A(g — r) ~ 0.1 redder than type II AGN, smaller than the 
binsize used. The colour distributions of spiral and elliptical 
AGN hosts differ only slightly; elliptical hosts are on av- 



erage only A(g — r) ~ 0.2 redder than spiral hosts. This 
small difference is common in AGN hosts, as they often lie 
in the "green valley" o f the colour-magnit ude relation dia- 
grams (see for instance iPadilla et al.|[2009h . 

We have thus constructed normal galaxy samples that 
reproduce the colour-luminosity-concentration properties of 
AGN populations (selected by AGN type and host morphol- 
ogy) which, unlike AGN of a given type, are expected to be 
randomly oriented with respect to the plane of the sky. We 
now characterise the intrinsic typical shapes and orienta- 
tions of the selected galaxy control samples. 

2.3 The intrinsic shapes of AGN galaxies in the 
SDSS 

We measure the intrinsic galaxy shapes in each AGN sample 
by studying the projected b/a distribution in the r-band for 
the control samples. 

Fig. [2] shows the distribution of minor to major pro- 
jected axis ratios (b/a from exponential or de Vaucouleurs 
fits corrected for the effects of the PSF) for spiral (top) and 
elliptical (bottom) type I (left panels) and type II (right 
panels) AGN galaxies (solid lines) and the corresponding 
control samples (dashed lines). All the distributions shown 
are Vinax- weighted. Errors are calculated using the jackknife 
technique. Both type II AGN ellipticals and spirals seem to 
be slightly more edge-on than their control samples, which 
are assumed to cover the full range of inclination angles. The 
shift in the b/a distribution peak towards lower b/a is about 
« 0.05 and « 0.15 for the spiral and elliptical type II popula- 
tions, respectively. There are considerably larger differences 
between the b/a distributions of the spiral type I AGN popu- 
lation and its control sample; the AGN population is skewed 
towards high b/a values, indicating inclination angles closer 
to face-on orientations. The control sample shows a peak in 
the distribution at b/a ~ 0.35, while the AGN population 
peaks at b/a ~ 0.65. This suggests that the source of ob- 
scuration in type II objects has a large scale height, and is 
preferentially aligned with the plane of the galaxy disc. 

At first glance, the larger inclination effect for type I 
AGN than that observed for type II AGN could be taken 
as a failure of the unified AGN model, since it naively sug- 
gests that the relation between the host galaxy and its AGN 
depends on the AGN type. However, a proper comparison 
between the two types of AGN requires a more detailed anal- 
ysis of relative inclination angles which we will perform in 
Section [3] another possibility that we will also explore is 
that the samples of different types of AGN are subject to 
different selection effects. 



2.3.1 Characterisation of the intrinsic shapes 

We use the b/a distributions to characterise the three- 
dimensional shapes of galaxies i n each control sample. A 
variety of papers (iRydenl |2004 IVincent fc Rvdenl 120051 . 
IPadilla fc Strauss!l2008l ) have shown that both elliptical and 
spiral galaxies have shapes consistent with oblate spheroids. 
In this work we model galaxies as triaxial ellipsoids of major 
axis A, middle axis B and minor axis C parametrised by two 
axis ratios, B/A and C/B. We follow the PS08 assumption 
that the distribution of 1 — C/B follows a Gaussian with 
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Figure 1. V ma x-weighted distributions of rgo/rso (left panels), Mr — 51ogio h (middle-left panels), g — r colour (middle-right), and 
redshift distributions (right panels) for AGN (solid lines) and control galaxy samples (dashed lines), divided according to AGN and 
morphological type. Absolute magnitudes and colours were corrected for internal extinction and reddening, respectively. 



mean 7 and standard deviation <r 7 , and the distribution of 
e = ln( l — B/A) is G aussian with mean /1 and dispersion a. 
Indeed lRvdenl (|2004l ) found that the distributions of intrinsic 
axis ratios are consistent with the Gaussian and log-normal 
shapes we are using. We use the PS08 model to fit the ob- 
served b/a projected distribution in order to determine these 
parameters for each galaxy sample. This approach also re- 
quires a model for internal extinction, as we described above. 
An important quantity is ip(8), the ratio of the number of 
galaxies seen at inclination 9 to the number expected with- 
out extinction, 



jrj^fE(M)f R (C) ( p(M)4>(C)W{C,M)dCdM 

mmcw(c, M)dMdc 



(i) 



where cj>(M) and (j)(C) are the uncorrected luminosity (here 
in the r-band) and colour (here g — r) functions, respectively, 
for a given sample. The function /e(M) = 4>e(M) / '4>{M) is 
the ratio between the number of observed extincted galaxies 
and intrinsic number of galaxies at a given luminosity, where 



the subindex E indicates that it refers to the extincted lu- 
minosity function. Similarly, /r(C) is the ratio between the 
underlying and reddened distributions of galaxy colours. W 
contains the correlation between g — r and M r , which we 
assume to be Gaussian with mean and dispersion as mea- 
sured directly from the data. The function ip(9) is thus com- 
pletely characterised by the dust-corrected luminosity and 
colour functions, and by the parameters of the dust and 
shape model, namely Eo, fi, a, 7 and rr 7 . A grid of param- 
eters is constructed, and for each grid point, pi, we obtain 
a model Af mo d e i(b/ a > {p}i) projected axis ratio distribution, 
which is compared to the observed distribution; 



x 2 aph)= E 

b/a bins 



N modcl {b/a, {p}j) - N{b/a) 

CTjack-knifc 



(2) 



The best-fitting parameters correspond to the minimum 
value of x 2 throughout the parameter grid. Following PS08, 
we assume Eo — for elliptical galaxies. 
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Figure 2. Distributions of projected axis ratio, 6/a, for selected spiral galaxies (left) and elliptical galaxies (right) hosting a type II 
(left panels) and type I (right panels) AGN (solid lines) and the corresponding control sample (matched to the distribution in the 
colour-luminosity-conccntration distributions that the AGN population follows; dashed lines). All the distributions have been weighted 
by 1/Vmax, corrected for extinction. Errors were calculated using the jackknifc technique. 



Table 2. Best-fitting model for each of the four control galaxy samples selected to mimic the colour, luminosity and concentration 
distributions of spirals and elliptical Type I and II AGN. The last column shows the likelihood associated with the parameter set. 



Sample Eq fi a 7 a 1 -Pmax 



Type I ellipticals 




0.0 


-0.9 


± 


0.5 


2.3 


± 


0.6 


0.45 


± 


0.03 


0.21 


± 


0.04 


0.86 


Type II ellipticals 




0.0 


-1.35 


± 


0.4 


1.7 


± 


0.5 


0.45 


± 


0.03 


0.23 


± 


0.04 


0.88 


Type I spirals 


0.3 


± 0.3 


-0.85 


± 


0.35 


1.7 


± 


0.2 


0.75 


± 


0.02 


0.07 


± 


0.03 


0.25 


Type II spirals 


0.3 


± 0.3 


-0.25 


± 


0.3 


2.2 


± 


0.2 


0.75 


± 


0.02 


0.04 


± 


0.02 


0.33 



Fig. [3] shows the best-fitting model compared to the ob- 
served distributions for each control galaxy sample. Sym- 
bols represent the observed SDSS galaxy b/a distribu- 
tions, and solid lines are the best-fitting models. The best 
fit parameters for each population are given in Table [2] 
The likelihood, P max , that corresponds to the minimum 
X 2 value quantifies whether our model is a good fit to 
the data, and is given for each case. The model axis dis- 
tributions follow the observed ones well. The values of /1 
and 7 for both elliptical galaxy populations are consistent 
with oblate sphe r oids, in agreement with previous analyses 
(e.g. lRvdenll2004l ; IVincent fc Rvdenll2005l ; iPadilla fe Strauss! 



|200S| ; iBattve fe Brownell2009l ). Moreover, control samples of 
both morphological types are characterised by shapes con- 
sistent (within the errors) with those of the full galaxy 
population as measured by PS 08 using the SDSS DR6 
ijAdelman-McCarthv et all fioOrj ). The measured distribu- 
tions, as well as the best fit values of Eo, 7 and ct 7 , are similar 
for the type I and II AGN spirals; fi and a show differences 
which could be used to argue against the unified model, 
but these are only marginally significant (slightly above one 
standard deviation; errorbars are given in Table These 
slight differences could be due to different selection effects 
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Figure 3. Distributions of projected axis ratios, (b/a), weighted 
by 1/Vmax for the control samples of galaxies selected to populate 
the same colour-luminosity-concentration space as our four classes 
of AGN. Symbols correspond to the control samples selected using 
the SDSS DR7 spectroscopic sample, and lines are the best-fitting 
models for the parameters shown in Table [2] 



acting on each AGN type, an issue that we will explore fur- 
ther in Section f3.ll 

In order to analyse possible parameter degeneracies we 
study the probability contours (calculated from \ 2 ) m the 
H — 7 and Eo — /i spaces. Fig. [4] shows n — 7 probability 
contours for the control samples of each AGN sample (as 
labelled). The white region surrounding the contours repre- 
sents probabilities < 0.01. The values of and 7 are very 
well constrained, particularly for spiral galaxies, which show 
a very well defined region of high likelihood, while for ellip- 
tical galaxies the contours tend to be slightly elongated and 
tilted, allowing a larger range of possible parameter values 
in /1. This may be due to the lower sensitivity of the method 
to large negative values of /1, which are associated with very 
round shapes. 

We also inspected the fi — a probability contours for 
all control samples (not shown) and find that these two pa- 
rameters are slightly more degenerate than the y, — 7 plane. 
The quantity a smooths the N mo< ici(b/a) distributions, and 
therefore, the relatively large uncertainty in this parameter 
simply indicates that it does not have an important effect 
on the shape of the modelled b/a axis ratio distributions. 

Fig. [5] shows the probability contours in the Eo — /j, 
space for the control samples of the type II (top panel) and 
type I (bottom panel) AGN spiral galaxy samples (ellipti- 
cals are assumed to have Eo = 0). The extinction is only 
marginally detected in both cases. Dust strongly determines 
the likelihood that a galaxy will enter the sample at a given 
viewing angle, which consequently affects the observed axis 
ratio distribution. Note that /1 is very insensitive to Eo, ex- 
cept for extreme values of /i (i.e. B/A <C 1). This means 
that when galaxies are intrinsically elongated (independent 
of their inclination angle), a large extinction is needed in 
order to broaden the projected b/a distribution and move 
its peak towards larger b/a. 




-1.5 



-2.0 { 



control sample 
Type II spirals 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
extinction 



0.00 

0.04 
0.07 
0.11 
0.14 
0.18 
0.21 




control sample 
Type I spirals 

-2.0 1 i i i i i i i i i i i i i i i i i i i i M 0.47 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
extinction 

Figure 5. Probability contours in the Eo- fJ. space, for the control 
samples of AGN spiral galaxy samples (as labelled) . The colour- 
scale is shown on the right of each panel. 



The modelled b/a distributions do not show significant 
changes within the errors calculated for each parameter, im- 
plying that the confidence regions are reasonable. Now that 
we have characterised the intrinsic shape distribution of each 
control sample, we use these measured shapes to estimate 
the viewing angles of the AGN populations as a whole, and 
for different AGN and morphological types separately. 



3 INCLINATIONS AND ORIENTATION 

ALIGNMENTS OF THE AGN POPULATION 

We aim to obtain the inclination angle distribution of the 
AGN host galaxy population to determine whether there 
is a preference for face- or edge-on orientations in a given 
AGN population. To do this, we first calculate the predicted 
likelihood distribution of the projected b/a as a function of 
the inclination angle (8), given our model for the intrinsic 
shapes of galaxies of each control sample (i.e. parameters 
listed in Table [2]). This is shown in Fig. [6] Here, the colour 
bar shows the likelihood normalized to give an integral of 
1 at each value of b/a. Elliptical galaxies are characterised 
by wide distributions in 9 at a given b/a as a result of their 
intrinsic round shapes, while spiral galaxies have much nar- 
rower distributions, supported by the value of the probabil- 
ity peak in each case. Note that there are two local maxima 
in the top-right and bottom-left of each panel. These repre- 
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sent a non-null probability of having an intrinsically round 
or elongated object which has an edge- or face-on inclina- 
tion, respectively. This probability is more significant in the 
case of elliptical galaxies. 

For each control galaxy sample we determine their dis- 
tribution of inclination angles using the distribution of mea- 
sured b/a values and the probability function of inclination 
angles vs. b/a shown in Fig. [6] These probability functions 
peak at higher values of inclination angle when b/a de- 
creases, and are narrower at lower b/a values in all the con- 
trol samples. We count the number of galaxies at each incli- 
nation angle, 6, weighting by W(9) = x Wgbi{0 \b/a), 
where Wg&i(8 \b/a) is the probability that a given galaxy 
has an inclination angle of 9 given its b/a. When we do 
this for the control samples we find that they are consistent 
within the errors with uniform distributions in cos(#), as ex- 
pected. To quantify the deviations from uniformity for the 
AGN sample, we study the ratios between the normalised 
distributions of cos(#) of the AGN population and the cor- 
responding control sample, /agn //control- 

Fig. [7] shows the distributions of /agn //control for spi- 
ral (solid) and elliptical (dashed) type II (top) and type I 
(bottom) AGN populations. The errorbars were calculated 
using the jackknife technique. Spiral type I AGNs have a 
clear tendency to be face-on; the tendency is weaker but 
still significant in the elliptical type I population. Type II 
ellipticals show the opposite tendency, with a distribution 



that has a maximum at 9 « 90°, while type II spirals show 
no clear skew in their inclinations. The weak signal in the 
latter cases is consistent with the similarity in the b/a distri- 
butions of the AGN and control samples (i.e. A b / a ~ 0.05, 
Fig. [2]), which can be interpreted in our model as differences 
in the polar viewing angle distributions. In general, the re- 
sults from Fig. [7] allow us to rule out random orientations 
(i.e. a flat distribution in /agn //control) with a confidence 
of 5 x 2 ~ 230 in the type I spiral population, S \ 2 ~ 3 in the 
type I elliptical population, S \ 2 ~ 15 in the type II elliptical 
population, and 5 \ 2 w 6 in the type II spiral population. 



There have been several previous attempts to study 
the b / a distributions of type I and II AGNs (e.g . iKeell 



1980, Lawrence & Elvis 




1982, 


Kirhakos fc Steinerl 


1990, 


McLeod & Riekc 1995, Naear & Wilson 1999, Riebv et all 



2006), finding that type I AGN usually appear to have axis 



ratios close to 1, in agreement with our findings. Our ap- 
proach is the first to include a consistent model that takes 
into account the intrinsic shapes of the AGN hosts and their 
intrinsic dust extinction. 

In the unified model the type of AGN one observes 
is just a matter of the inclination of th e obscuring torus 
to the line-of-sight (e.g. lAntonuccilll993l ). Within this sce- 
nario, skewness in the distribution of the inclination angles 
of host AGN galaxies, like the ones we have found here, 
would reflect a preferential alignment between the large 
structures of the galaxy (kpc scales) and their central re- 
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Figure 6. Probability density contours in the axis ratio, b/a, vs. inclination angle plane for each galaxy population selected to mimic 
the luminosity, colour and concentration properties of the AGN population indicated in each panel. These results take into account the 
effects of dust extinction. The scale is shown on the colour bar at the right of each plot. Probabilities were normalized to give an integral 
of 1 in horizontal cuts of b/a. 



gions, as one would expect in the coherent scenario (see 
LPC09). However, the galactic disc itself may contribute 
to the absorption of the emission f rom the central engine , 
and particularly from the BLR (e.g. iMaiolino fc RiekeJ | l995 |; 
Lutz et al.l I200&1 ; IMaiolino et~atl 120031 ; ISatvapal et al-lhool 



Goulding fc Alexanderl l2009 . 2010), giving rise to skewness 



in the inclination angle distribution of AGN hosts. This is 
why it is important to explore the nature of the observed 
tendencies towards face- or edge-on of Fig [7] 

In order to gain insight into the nature of the inclina- 
tion angle distribution of the AGN samples, we study the 
dependence of the frequency shown in Fig[7]on properties of 
the [OIII] emission line. We also calculate the expected dis- 
tributions of b/a for spirals if the galactic disc is responsible 
for the absorption of BLR emission in type II AGN. 



3.1 Possible completeness problems in the AGN 
samples 

Line emission from star formation activity in galaxies could 
represent an important bias in the optical classification of 
AG N from their emission-line p roper ties. Several authors 
fe.glGoulding fc Alexanderll2009l . 2010; I Juneaufc 7 Dickinson 



l201lh have shown that a significant population of Seyfert lis 
hosted by edge-on galaxies are misclassified due to this effect 
at low and intermediate redshifts. We expect this effect to be 



more important at low emission line luminosities since in this 
regime, the star-formation contribution to the emission-line 
flux could overwhelm the AGN light. Note that faint, low- 
luminosity emission lines are also difficult to detect given our 
high S/N limits. In this section we explore the sensitivity of 
the inclination angle distribution to the EW of the lines, 
with these effects included. 

We divide the four samples of AGN galaxies in subsam- 
ples of low and high rest-frame EW of the [OIII] emission 
lin<Q (divided at the sample median [OIII] EW of « 3 A). 
For each subsample we construct its corresponding control 
sample as in §2. We find that these control samples are char- 
acterised by the same concentrations and g — r colours as 
those shown in Fig. [T] However, the typical luminosities of 
the high and low [OIII] EW samples are brighter and fainter 
by ~ 0.3 magnitudes than the distributions of Fig.[T] respec- 
tively. Th is is consistent with the observed BH-bulge rela- 
tions fe.g.lFerrarese fc Merritt|[200ol ; iMarconi fc Huntll2003l ; 
lHaring fc R ix 2004), where more massive BHs (which drive 
brighter AGN, in this case linked with higher [OIII] EW) are 



2 We remind the reader that the [OIII] emission line is generated 
in the narrow line region of AGNs, and therefore always appears 
narrow independent of the AGN type. This enables us to do a fair 
comparison between type I and II AGN subsamples. 
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Figure 8. Ratios between normalised frequencies of cos(0) of the AGN population and the corresponding control samples, for the spiral 
and elliptical type II (top) and type I (bottom) AGN populations, separated according to their [OIII] equivalent width (low EW in solid 
lines and high EW in dashed lines). 



Table 3. Weighted fraction of edge-on galaxies in the type I 
and type II AGN spiral galaxies samples compared with the cor- 
responding control sample, and in subsamples of high and low 
[OIII] equivalent width (see Fig- [Sj • Edge-on galaxies are defined 
as those having an inclination angle of cos(6) ^ 2(1 — 7) (i.e. 
corresponding to a maximum deviation from the line of sight of 
1 — 7 =< C/B >, see Table [2{. Errors indicate the standard 
deviation of the mean. 



Sample 


AGN sample 


control sample 


All Type II AGN spirals 


0.58 ±0.05 


0.56 ±0.05 


Low [Olll] EW 


0.59 ± 0.07 


0.63 ±0.06 


High [OIII] EW 


0.70 ±0.03 


0.62 ±0.04 


All Type I AGN spirals 


0.28 ± 0.06 


0.44 ± 0.05 


Low [OIII] EW 


0.36 ± 0.06 


0.56 ±0.07 


High [OIII] EW 


0.43 ± 0.07 


0.62 ±0.05 



hosted by more massive galaxies. Despite this difference in 
the intrinsic luminosity of the hosts, the intrinsic shape pa- 
rameters derived from the control samples in each subsample 
of low and high [OIII] EW agree with each other within the 
uncertainties, and with the parameters found for each AGN 
sample overall (Table [2}. Interestingly, the marginal differ- 
ences seen between the intrinsic shapes of type I and II AGN 
spirals in Table [2] disappear when comparing the high [OIII] 
EW subsamples. 



The resulting cos(#) distributions of the subsample of 
AGN galaxies of high and low [OIII] EW are shown in Fig.! 
The low/high [OIII] EW subsamples are in good agreement 
with one another within the same AGN sample, except for 
type II spiral galaxies, where the high [OIII] EW subsample 
shows a strong edge-on tendency. 

Since the control samples have underlying shapes in 
agreement with each other, the differences seen between the 
low and high [OIII] EW subsamples of the same AGN sam- 
ple in Fig. [H] are not due to intrinsic differences in the hosts. 
This suggests that the low EW [OIII] subsample of type II 
spirals is subject to larger systematic biases than its high 
EW counterpart, and the unified model does provide a good 
explanation of the results for the latter. 

In order to quantify the differences between the low and 
high [OIII] EW subsamples of type II AGN spirals, we calcu- 
late the VMAx-weighted fractions of edge-on galaxies (those 
with cos(0) < 2(1 - 7) = 2 < C/B >, i.e. 9 > 60°) in each 
subsample. These fractions are shown in Table [3] For ref- 
erence, we also show the edge-on galaxy fraction for type I 
AGN spirals, and for low and high [OIII] EW subsamples of 
the type I spirals. This fraction for the high [OIII] EW sub- 
sample is 20% higher than in the low [OIII] EW subsample 
in type II spirals, and is larger than the errors of the means. 
This suggests we are missing m 20% of edge-on type II spi- 
rals with low [OIII] EW. This number is a lower limit, since 
it assumes that the high [OIII] EW subsample is complete 
and that we are just losing edge-on galaxies. This crude es- 
timate is consistent within a factor of 2 with the estimate 
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Figure 7. Ratios between normalised frequencies of cos(6) of 
the AGN population and their corresponding control samples, 
/agn //control > f° r the spiral (solid lines) and elliptical (dashed 
lines) type II (top) and type I (bottom) AGN populations. The 
distributions are 1/Vmax-weighted. Errorbars were calculated us- 
ing the jackknife technique. 



by I Juneau &: Di ckinson 20j3] of optically misclassified SDSS 
DR4 AGN identified by their X-ray emission. They found 
that 40% of X-ray-bright edge-on galaxies are classified as 
star forming or composite star forming/AGN galaxies. 

In the case of type I AGN spirals, the high [OIII] EW 
subsample shows a larger fraction of edge-on galaxies than 
the low [OIII] EW subsample, possibly indicating that the 
bias mentioned above is also affecting the type I AGN sam- 
ple. However, when compared to the control sample, the dif- 
ference is only 8%, which is smaller than the errors. Thus, 
the two subsamples are consistent within the errors, as ex- 
pected from the distributions of Fig. [3] We conclude that 
obscuration from the galactic disc only has a significant ef- 
fect on our ability to classify type II spirals. 



3.2 Could the broad line absorption come from 
the galactic disc? 

The obscuration preventing the direct observation of the ac- 
cretion disc in type II AGN galaxies could be produced by 
a local torus or the galactic disc. Even though polarimetry 
results support the unified AGN model which invokes the 
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Figure 9. Expected distributions of b/a when the absorption 
of the broad emission lines comes from the disc of the galaxies 
instead of a nuclear torus, normalised to the number of galaxies 
in the samples of AGN galaxies (dashed lines). The distribution 
of projected axis ratios, b/a, for the type II (top panel) and type 
I (bottom panel) AGN spirals are shown as solid lines, and the 
corresponding control samples as dotted lines. For reference, we 
also show as a dot-dashed line the composite b/a distributions for 
the case where the galactic disc is half as thick (i.e. smaller C/B). 



local torus, detailed analysis of nearby galaxies have also re- 
vealed absorpt ion of broad line regions by t he galactic disc 
m some cases (|Goulding fc Alexanderll2009l ; iGoulding et al.l 
|2010| ). 

The results of the previous sections do not distinguish 
between these two possibilities. In an attempt to address this 
question, we test whether the observed b/a distributions of 
type I and II spiral AGN can be reproduced by absorption 
by dust in the discs of their host galaxie^. To do this we use 
the heuristic dust model adopted by PS08, which depends 
on the parameter 7, the mean height of the galaxy ellipsoid. 
We assume that the dust lanes produce a covering fraction 
which depends on the luminosity of the galaxy and that 
varies with the inclination angle 9 of the galaxy as 



3 It has been pointed out that elliptical galaxies often show cir- 
cumnuclear discs w hich can also obscure the central engine (e.g. 
iKawata et al.ll2007l) . However, our method is not able to probe 
this possibility. 
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/cov(#) — 



fo x 10°- 4(7 - cose) ,if COS* 
/o, if cos0 < 7 



> 7 



(3) 



of the solid angle covered by the torus corresponds to the 
weighted ratio of type II AGN to the total number of AGN, 



where fo is the edge-on covering fraction. Comparing with 
§2, the io°- 4 <t- cos e ' factor reflects the fact that we are con- 
sidering extinction in luminosity rather than magnitudes. 
We adopt the extreme case of fo = 1. We produce distri- 
butions of b/a for galaxies with inclination angles in bins 
of cos(#) weighted by the corresponding covering fraction 
/cov, add them together, and compare the result with the 
observed b/a distribution of type I and II AGN spirals. 

This is shown in Fig. where the observed b/a dis- 
tributions for the AGN samples (see Fig. [2]) are shown as 
solid lines, and the composite distributions for the case of 
pure disc absorption (fo = 1) are shown as dashed lines. 
For reference, we also show the composite distributions in 
the case where we consider galactic discs that are narrower, 
i.e. in C/B, by a factor of two. The latter choice is moti- 
vated by the fact that the PS08 model attempts to fit the 
bulge and the galactic disc as a single photometric structure 
instead of distinguishing between them, which could spuri- 
ously broaden the galactic disc. Absorption by the galactic 
disc alone goes in the right direction, but does not fully re- 
produce the distributions of either type I or II AGN even 
assuming fo = 1; the resulting distributions show excesses 
of face-on and edge-on galaxies in the two cases, respec- 
tively. By considering narrower galactic discs we are able to 
broaden the b/a distributions for type I spiral AGN slightly. 
However, the output distribution still has a clear excess of 
face-on galaxies. Moreover, / cov is not strongly dependent 
on 8 (i.e. a large population of type II AGNs could come 
from fairly face-on galaxies), meaning that we are forcing 
the composite b/a distribution to be as broad as possible. 
Thus a more physical approach would give even narrower 
distributions. In addition, the pure galactic disc absorption 
scenario adopted here produces five times more type II than 
type I AGN at a given lumin osity, in contradic tion with the 
observed ratio of 0.25 - 0.^3 i|Hao et al.ll2005bl ; see also §4). 

Our results indicate that even if part of the broad line 
absorption takes place in galactic discs, additional sources 
of absorption such as a torus, are still needed. We will bear 
this in mind when analysing the relative orientations of torus 
and galactic disc angular momenta in the following section. 



4 ALIGNMENTS BETWEEN THE GALAXY 
DISC AND AN OBSCURING TORUS 

We use the measured inclination angle distribution (cf. 
Fig. to analyse alignments between the galaxy disc and 
the broad torus in spiral galaxies within the framework of 
the AGN unified model. 

Within this paradigm, the obscuring torus must be 
broad in order to explain the wide inclination angle range 
over which type II AGN are found, and the ratio b etween the 
number of type II and type I AGN galaxies ( e.g. lHao et all 
l2005bl ; iTreister et afll200Sl : iReves et alj|2008h . We quantify 
the typical torus width by considering that the fraction 



r i = 1 . 



Yl Wagntti 
J2 Wagnti ' 



(4) 



where Wagn = Kiii- This calculation shows that the torus 
has a typical azimuthal height of ~ 40°, which corresponds 
to a ViviAx-weighted ratio of type I to II of 3:1. By measur- 
ing the [OIII] luminosity function of type I and II AGNs, 
Hao et al. (2005b) show that this ratio depends on the AGN 
luminosity. They find ratios close to unity for low luminosi- 
ties, ~ 2 : 1 for 10 5 - 8 < L([OIII]) < 1O 65 L , and ~ 4 : 1 
for L([OIII]) > 10 6 ' 5 Lq, broadly consistent with our esti- 
mate, considering that the weighted mean luminosity of our 
sample is L([OIII]) « 10 6 4 L Q . 

We use the output galaxies from the semi-analytic 
model of LPC09 to construct samples of modeled type I and 
type II AGN spirals and ellipticals using the typical torus 
width inferred above. We choose a random observer's line of 
sight and estimate the projected angle subtended by the an- 
gular momentum of the accretion disc, # a cc- We assume per- 
fect alignment between the obscuring torus and the accretion 
disci We run a Monte-Carlo simulation and extract torus 
heights, ht, so that the mean corresponds to tan(40°) w 0.83. 
We consider those galaxies with 6> acc + h t > tv/2 as type II 
AGN, and the rest as type I AGN. Given each galaxy angu- 
lar momentum direction, we estimate the observed projected 
inclination angles. 

The left and middle panels of Fig. [TO] show the result- 
ing normalised frequency of the inclination angle of AGN 
hosts for the LPC09 semi-analytic model in the coherent 
and chaotic scenarios, and compare this with the observed 
distributions from Fig. [7] For reference, we also show the 
inclination angle distributions for the high [OIII] EW sub- 
samples from Fig. [H] since they are more likely to be com- 
plete samples (see §3.1). The right panel of Fig. [10] shows the 
theoretical expectations for the angle between the angular 
momentum of the galaxy and the accretion disc in the co- 
herent and chaotic scenarios. In the chaotic model, the two 
angular momenta are unrelated and the distribution is flat. 
In the coherent model they are highly aligned due to the 
assumption that the accretion flow to the BH is perfectly 
aligned with the large-scale source of the gas. The coher- 
ent scenario reproduces the observed distributions of type I 
AGN quite closely, but it fails for type II AGN (particularly 
for type II AGN spirals) . The discrepancies decrease slightly 
when comparing to the high [OIII] EW subsamples. 

As was mentioned above, obscuration by the galactic 
disc may also affect the determination of the AGN type in 
spiral galaxies (see §3). We consider this possibility by also 
including the heuristic model for the obscuration of the AGN 
emission by the galactic disc of §3.2 in the predictions of the 
LCP09 model in the coherent and chaotic scenarios. In this 
case, highly inclined spiral galaxies have higher probability 
of being seen as type II AGN. A fraction / CO v(#) of spiral 
galaxies at inclination 6 will be classified as type II AGN. 
Fig. [11] shows the distribution of cos(#) for spiral galaxies 



4 Note that these ratios are VMAX _we ighted and thus refer to 
weighted number of detections. 



5 This is expected at least to a high degree in models where 
the absorption of BLR occurs i n radiation-driven w inds from the 
accretion disk (see for instance iMurrav et al.|[l995h . 
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Figure 10. Left and middle panels: Distribution of the inclination angle of AGN host galaxies from Fig. [7] (in solid lines; AGN types 
are as labelled) compared to the predictions of the semi-analytic model LPC09 in the coherent (dot-dashed lines) and chaotic (dotted 
lines) scenarios. For reference, we also show the cos(#) distributions for high [OIII] EW subsamples of Fig. [8] as dashed lines. Right 
panels: Predicted distributions of the angle between the angular momentum of the galaxy and the acc retion disc for the LPC09 in the 
coherent (solid line) and chaotic (dashed line) scenarios for elliptical and spiral galaxies ( classified as in lLagos et al.ll20o3 to account for 
the observed morphology type fractions as a function of stellar mass from lConse licc 2006). Error bars were calculated using the jackknifc 
technique. 



from Fig. 1101 together with the chaotic model and the co- 
herent model when the obscuration from the disc is included. 
The tendency of alignment observed in the LPC09 model in 
the coherent scenario becomes weaker for the type II spirals, 
and stronger in type I spirals. A tendency towards edge- and 
face-on orientations appears for type II and I AGN, respec- 
tively, in the chaotic model. This tendency is very similar 
to that observed: a very weak edge-on tendency for type II 
spirals, and a much stronger tendency towards face-on ori- 
entations for type I spirals. This suggests that the observed 
stronger signal for a preferred inclination of type I com- 
pared to type II hosts may be driven by obscuration from 
the galactic disc. However, the model of galactic disc obscu- 
ration we have used is maximal in the sense that it assumes 
that edge-on AGN galaxies are observed as type II AGN. 
More importantly, the adopted disc obscuration cannot re- 
produce the distribution of projected shapes of AGN hosts 
(see Fig. [9]) or explain the tendencies observed in the ellipti- 
cal population. Thus, disc obscuration cannot fully explain 
the type I and II dichotomy reported here. 

On the other hand, the coherent model predicts a higher 
skewing towards face- or edge-on orientations of AGN hosts 
than we see in the observations, indicating that the degree 
of alignment between the galaxy and the torus (or accretion 
disc), in the observational sample is much weaker than in 
the idealized coherent scenario. 

We summarize the comparison between the observed 
cos(#) distribution of each AGN sample and the predictions 



of the LPC09 model with and without disc obscuration in 
Table 3] As can be seen from the 5x 2 between observations 
and models, neither, the chaotic nor the coherent model, 
is able to closely reproduce the observed distributions. The 
only exception is the high [OIII] EW subsample of type I 
AGN ellipticals, which is consistent with random orienta- 
tions. This suggests that in general, the coherent and the 
chaotic scenarios represent idealized conditions; i.e. some 
loss of the angular momentum direction is expected during 
the infall, but this cannot be complete. 

Our results indicate that, within the AGN unified model 
framework, there is some weak coherence between the torus 
and the host galaxy angular momenta, and that the purely 
chaotic scenario is ruled out. 



5 DISCUSSION 

The results shown in the last two sections point to a scenario 
in which the AGN components, i.e. torus and accretion disc, 
are aligned to a significant degree with the angular momen- 
tum of the host galaxy. 

Theoretical suggestions by LPC09 indicate that, under 
this condition, high BH spin values would be commonplace 
(cf. Section 1). This preference for alignments could be ex- 
plained by short timescales chara cterising the gas i n fall to 
the central parts of the galaxy. iBogdanovic et all (|2007l ) 
show that this is likely to happen in gas-rich mergers (see 
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Table 4. 8\ 2 = X 2 /d.o.f. from the comparison between the inferred distribution of inclination angles of AGN hosts in the SDSS and the 
predictions of the LPC09 model in three variants: (i) the chaotic scenario, (ii) the coherent scenario (see Fig. HOP , and (iii) the coherent 
and chaotic scenarios considering obscuration from the galactic disc (only applied to spiral galaxies; see Fig II IB . 
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Figure 11. Distribution of the inclination angle of AGN host spi- 
ral galaxies from Fig. [7] (in solid lines; AGN types are as labelled) 
and for the high [OIII] EW subsamples of Fig. [8] (dashed lines), 
compared to the predictions of the semi-analytic model LPC09 
in the coherent (dot-dashed lines) and chaotic (dotted lines) sce- 
narios, when the effect of galactic disc obscuration, by using the 
heuristic model of 53.2, is included. 



also iBarausse fc Rezzollalliurol ). where the BHs taking part 
in the galaxy merger will acquire 1%-10% of their mass in a 
time which is short compared to the time neede d for the BHs 
to sp iral in towards the center 5 x 10 7 yr, Esca la et al.l 
|2004| . 2005), or t he time for a starburst to deplete the supply 
of gas (~ 10 8 yr, Larson 1987). A potential problem with the 
preservation of the angular momentum direction in inflow- 
ing gas assumed in the coherent model comes from the dif- 
ficulty in transp orting gas through th e corotatio n resonance 
radii of galaxies i|Zhang fc Butall2007l ). However. lHaan et al.l 



(2009) found evidence for gas inflow b eyond this radius in 
seven galaxies in the sample studied bv lGarcia-Burillo et al.l 
(2003), which continues to their resolution limit, lOOpc away 
from the nuclei. 

Direct comparis on with detailed simulations of the 
feeding of BHs (e.g. iHobbs et alj l20ld : iPower et al.ll2010l ; 
iHopkins fc Quataertl [20I0I ) represents an attractive tool to 
test accretion scenarios and compare their expectations with 
the orientation results presented in this work. Nonetheless, 
BH feeding is indeed a challenging issue; understanding 
it depends strongly on modeling the local environmental 
effects such as star formation act ivity, supernovae explo- 
sions, local winds, and so on (e.g. ICiotti fc Ostrikerll2007l ; 
iKing et al.ll200Sl ; iNavakshin fc Powerll2010l) . ~ 



6 SUMMARY AND CONCLUSIONS 

We study the intrinsic shapes of SDSS DR7 AGN host galax- 
ies and their inclination angle distributions, under the as- 
sumption that AGN host galaxies have the same underlying 
shapes as normal galaxies matched in g — r colour, r— band 
luminosity, and concentration, rso/rgo. We examine ellipti- 
cal and spiral type I and II AGN galaxies separately, and 
ch aracterise the 3-D s hapes of these samples using the model 
of iPadilla fc Strauss! (|2008l ). With these models we infer the 
distributions of inclination angles of AGN hosts. The main 
results are as follows. 

(i) The structural parameters of the AGN control sam- 
ples (Table 1) are consistent with the fu ll SDSS DR6 spi- 
ral an d elliptical population reported bv lPadilla fc Strauss! 
(2008), and correspond to oblate spheroids with typical 

7 =< 1 — C/B > of « 0.40 ± 0.24 for elliptical galaxies, and 
« 0.77 ± 0.08 for spiral galaxies. The unified AGN model 
states that type I and II AGN are similar objects seen with 
different orientations. Consistent with this, we find that the 
intrinsic shapes of type I and II AGN hosts of a given mor- 
phological type are comparable. 

(ii) Using the intrinsic shape parameters of control sam- 
ples, we find that type I AGN galaxies have a strong ten- 
dency to be face-on, while type II AGN galaxies have only 
a slight tendency to be edge-on. In particular, type II spiral 
galaxies show a more uniform inclination angle distribution 
consistent with random orientations; only galaxies in the 
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subsample of high [OIII] EW show clear edge-on orienta- 
tions, indicating a possible selection effect acting on the low 
[OIII] EW spiral AGN. In the case of type I galaxies we are 
able to rule out random orientations with a confidence of 



sx 2 



3 and 8\ ~ 230 in the elliptical and spiral pop- 



ulations, respectively. In the case of the type II elliptical 
and spiral galaxies, we rule out random orientations with a 
confidence of 8\ 2 ~ 15 and S\ 2 ~ 6, respectively. 

(iii) We use the estimated three-dimensional shapes of 
spiral galaxies to test whether dust extinction in the galactic 
discs could be responsible for the absorption of broad AGN 
lines in type II AGN galaxies. We find that the resulting 
predicted b/a distributions are not compatible with those 
observed for either type I or II AGNs. This indicates that 
galactic discs cannot be the only source of broad line absorp- 
tion, and supports the existence of broad tori surrounding 
the active nuclei. 

(iv) Using the weighted frequency of type I to type II 
AGN, we found that the torus producing the absorption of 
broad lines has a typical azimuthal height of ~ 40 degrees. 

(v) We compare the observed inclination angle distri- 
butions with the theoretical predictions of the LPC09 semi- 
analytic model for the coherent scenario under the AGN 
unified model (e.g. iFanaroff fc Rilevl Il974l ; iBarthell 1 19891 ; 
iMadau et al1ll99-tlGunn fc Shankslll999l ). We find that the 
model overpredicts the observed skewing towards face- or 
edge-on orientations of AGN hosts, but sucessfully repro- 
duces the differences in the skewing between ellipticals and 
spirals. On the other hand, the chaotic scenario predicts ten- 
dencies towards face- and edge-on orientations of the type 
I and II AGN spirals, respectively, only if absorption from 
the galactic disc is assumed. However, in the case of ellipti- 
cal galaxies the chaotic scenario fails unless large amounts of 
dust aligned with the major axis of the galaxy are assumed. 
Our results suggest that some, but not all, of the direction of 
the angular momentum of the material falling to the nucleus 
is coherent. 

These results put important constraints on the physical 
processes involved in the gas inflow from the outer parts of 
the galaxy to the central engine. If the coherence of the ac- 
cretion flow is frequent and large enough, high spin values 
(a ~ 1) would be commonplace at the massive end of the 
BH population (|Lagos et al. 2009), regardless of other phys- 
ical processes as warps (IKing e t al. 2005) and fragmentation 
inside the accretion disc (|King et al.ll200Sl ). However, a sce- 
nario where just a small fraction of the gas flow keeps its 
original angular momentum direction has not been studied 
theoretically yet, and its consequences on BH growth and 
its spin remain to be explored. 

Kinemat i c stud ies of individual AGN galaxies, as in 
iDumas et all 1120071) . can help to further explore the align- 
ments between gas infall and the whole galaxy. It should 
also be borne in mind that our approach of fitting the in- 
trinsic shapes of galaxies from the photometry in SDSS has 
its limitations; a possible im provement could com e from us- 
ing the Galaxy Zoo Project (|Lintott et alj|2008h morpholo- 
gies, which can help to obtain more precise intrinsic shapes 
and, therefore, better and more reliable distributions of in- 
clination angles for AGNs of different types. A different and 
promising approach based on the study of the orientations 
of host galaxies of radio sources, selected as to indicate rel- 
ativistic jets pointing close to the line of sight, would give 



insight on the alignment between the galaxy and the BH 
spin (Lagos et al. in prep.). 

Looking furthe r to the future, the LISA survey 
l|johann et al.| [2008) along with new high-resolution X-ray 
spectroscopy, will provide unequivocal information on the 
typical spins of the BH population and the orientations of 
jets. This will be complemented with detailed kinematic 
studies of the gas in large samples of galaxies (revealing 
the gas inflows) th at the Square Kilometer Array (SKA; 
ISchilizzi et al.| [2008) will achieve, helping to solve these still 
open questions. 
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